Anomalous Hall Effect in Disordered Multi-band Metals 
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We present a microscopic theory of the anomalous Hall effect in metallic multi-band ferromagnets, 
which accounts for all scattering- independent contributions, i.e., both the intrinsic and the so-called 
side jump. For a model of Gaussian disorder, the anomalous Hall effect is expressed solely in 
terms of the electronic band structure of the host material. Our theory handles systematically the 
interband-scattering coherence effects. We demonstrate the method in the two-dimensional Rashba 
and three-dimensional ferromagnetic (III,Mn)V semiconductor models. Our formalism is directly 
amenable to ab initio treatments for a wide range of ferromagnetic metals. 
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Introduction. — While the anomalous Hall effect (AHE) 
has attracted generous attention from the physics com- 
munity starting with the seminal work by Karplus and 
Luttinger its full theoretical understanding remains 
incomplete |2j. Contributions to the AHE in ferromag- 
netic metals can be separated according to the depen- 
dence on the quasiparticle transport life time r, e.g., 



o-ah(t) = ay H + a H cj xx (T) + 
scattering- independent contribution, and <Xho~xx{t), usu- 
ally termed skew scattering, is linear in r in the Drude 
limit (i.e., nop 1, where fnop is the Fermi energy). 
The scattering-independent term o~^h ' 1S usually further 
separated into the intrinsic contribution (IC), o~™\i, and 
the side-jump contribution (SJC) aj 



where cr^ H is the 



= rr (0) 
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c™ H is defined as the extrapolation of the ac Hall conduc- 
tivity to zero frequency in a clean system, with the limit 
t -1 — > taken before u — > 0|. The IC has been shown 
to be linked to the Berry phase of the spin-orbit (SO) 
coupled Bloch electrons [3|. It is the most directly calcu- 
lable contribution to the AHE in ferromagnetic semicon- 
ductors, transition metals and complex oxides 0, @|- 

A wide range of strongly SO coupled ferromagnetic 
metals exhibit scattering-independent o~^ H in the a ah 
signal with a sizable deviation from the calculated <7™j^ 
USUI, which implies substantial SJC. Although an ex- 
perimental separation of IC and SJC has been suggested 
by studying an interplay of different kinds of disorder 
(e.g., phonons and impurities) at finite temperatures @, 
comparison with the theoretical expectation for cr^ H has 



is modified by the effects of disorder at low frequencies, 
while at intermediate frequencies, r _1 < u> < A, the IC 
should be recovered as interband coherences caused by 
disorder scattering do not build up Q. 

In this Letter, we calculate the AHE in metallic non- 
interacting multi-band systems in the presence of delta- 
correlated Gaussian disorder, expressing the final result 
solely through the Bloch wave functions, similar to the 
theory of the intrinsic AHE There is no skew- 

scattering contribution (oc o~ xx ) within such model of dis- 
order, and we assume nondegenerate bands. The main 
results of this Letter for the IC and SJC are given in 
Eqs. ©-O requiring only the material's electronic band 
structure as the input. These equations should apply to 
a wide range of metallic materials exhibiting scattering- 
independent AHE The present theory has been tested 
on the two-dimensional (2D) Rashba Hamiltonian repro- 
ducing known results [3, |9§. Furthermore, the SJC is 
found to dominate the AHE in a model of metallic ferro- 
magnetic (III,Mn)V semiconductors. 

The Berry phase of Bloch states has a significant ef- 
fect on transport properties, particularly on the AHE. 
The origin of this lies in the anomalous velocity pro- 
portional to the external electric field that modifies the 
group velocity i.e., hi — c^e^k) — eE x B v (\t), 
where e„(k) is the band energy, E external electric field, 
B v (k) = idk x (u^d^lurj) Berry curvature, and e par- 
ticle charge (e < for electrons). Modifications to the 
motion of electrons (holes) in the 77th band are defined 
solely in terms of the periodic part of the Bloch wave 



been hampered by the lack of a simple rigorous formalism functions |u,(k)). Below, we will show that this is no 
that would allow a reliable calculation of a s j H in complex longer the case due to band mixing, in the presence of an 

even infinitesimally small disorder. 

IC and SJC from the band structure. — Consider a gen- 
eral multi-band noninteracting system, in the kp expan- 
sion, up to some order in k about an extremum point in 
the Brillouin zone. Below, we will consider Luttinger and 
Rashba Hamiltonians as particular realizations of such 
expansions. In the position representation, our iV-band 



multi-band systems. It is thus desirable to develop a gen- 
eral procedure for calculating all scattering-independent 
contributions to allow for a systematic comparison with 
experiments and engineering of materials with necessary 
AHE properties. It should be possible to identify the SJC 
by ac measurements in the clean limit (i.e., t _1 < A, 
the characteristic SO band-energy splitting): The AHE 



2 



projected Hamiltonian is expressed via "envelope fields": of TV 2 linear equations with N equations 



Ho 



E 



drttt(r)[-ffo(-iV r )W*;, 



(27T) 
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Here, ^f v is the "envelope field" of the 77th band, with 
index rj running from 1 to N. We suppose that all in- 
formation about our system, such as SO interaction or 
exchange field, is contained in the matrix structure of 
Ho(k), where k corresponds to — iV r . The phenomeno- 
logical exchange field is introduced within the framework 
of a mean-field description. In addition to the band 
Hamiltonian, we include a scalar delta-correlated Gaus- 
sian disorder V{r) with (V(r)V(r')) = h 2 VS(r - r'). 

In the absence of disorder, the anomalous velocity men- 
tioned above leads to the intrinsic spin Hall conductance 



int e \ ~> 



d n k duj 
(2tt)™ 2?r 



npilA^Akj — AkjAktlrinAf, , 

. (2) 

where n = 2 or 3 is the number of spatial dimensions, 
tif(w) is the Fermi distribution function, and A«j(k, uS) 
is the spectral function of the 77th band. The anoma- 
lous transport is governed by the Berry-connection ma- 
trix Ak = iU'dkU, where U is a k-dependent unitary 
matrix transforming the k • p Hamiltonian into a diago- 
nal band-energy matrix e(k) = W ' H {k)U . 

In this Letter, we show that the SJC contains two 
terms expressed via the electronic band structure as 



N 

E 

x Tr 



1 



(27T) n |9ft k £r,| [TcJtjtj 

{(tfcUUAktSnU* - USr,A ki %S v U^ P 3 
+ (d hk] e v )S n A kt (i - §„)%} + c.c. , (3) 



N 



E ' 



57=1 

x Tr 



(2tt)™ 2|5 fik £ J7 | [%} m 



+ [ic] m US v U^P t UC v U^P 3 

+ 2(d hk ^)US n %C ri U^P 3 }+c.c. 



(4) 



Here, [5^]^ = Sijdi n (Sij is the Kronecker delta symbol) 
and [C n ]ij = (e v — £i)~ 1 5ij for i 7^ 77 and zero otherwise. 

•m - wit I r£, aoo , 



V/ 



(27T)™- 1 2|0 ftk e„| 



(5) 

and the matrix P corresponding to a subset of vertex 
corrections denoted by T in Fig. Q] is defined by a total 



for each 77. In the above equations, stands for 

the integration over the Fermi surface of the 77th band. 
The SJCs in Eqs. ^ and (gj) are distinct from the di- 
agrammatic point of view as will be clear below. The 
mechanism of the former SJC relies on the effects related 
to the Berry curvature hence the dependence on Ak,. 

Derivation. — In various theories of the AHE in multi- 
band systems, it is common to express the conductivity 
via the Green's functions (GF's) calculated in equilib- 
rium Such description requires as input information 
about both the disorder and band structure. By taking 
advantage of the band eigenstate representation, we will 
express our results only via the band structure. To fulfill 
this, we first express all GF's via their diagonal parts as 



U d ^nd 



i,R(A) 



G 



= G 



RCA) 



^R(A)^R(A)AR{A) 



(7) 



Here and henceforth, the eigenstate representation is de- 
noted by index c, E 



R(A) 



W^ A) U 



,R(A) 



,R(A) 
J nd 



is the self-energy matrix separated into the diagonal and 
off-diagonal parts, G^q A ^ = h(huj — H Q — Tjcq^)^ 1 is 
the retarded (advanced) GF in equilibrium and G^ A ^ = 
h(hjj — H — T,^^)^ 1 is the corresponding diagonal GF. 
The imaginary part of G^ A ^ proportional to the spectral 
function A = i(G^ — G A ) will be integrated out reducing 
the problem to Fermi-surface integration. It is crucial to 
keep off-diagonal matrices E^j in Eq. Q up to the nec- 
essary order (the first order for the leading-order AHE) 
since they contain information about the interband co- 
herences that play an important role in the AHE. 

Our starting point is the expressions for the current 
densities derived within the Kubo-Streda formalism (Io| 
by summing all noncrossed diagrams. These expressions 
are also obtained in [§| using the Keldysh formalism: 



Ji = - T 



d n k du 



EcLn F Tr 



VG^pGifii 



h J (27t)" 2tt 

(6* q vG A q - (G* q vG* q + G^vG")/: 



(8) 



d n k duj. 



h J (27t)™ 2tt 

R 



Era F Tr 



(Gf^cLG^ - d u G? q vG«) m/2 



where the vector- valued matrices p{u>) and p^(w) satisfy 
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the following equations (v = dnuHa)- 



P = 



VG*pG, 



G e qVG eq 



d n k 

(2tt)" 
d n k 
(2ir) n 



VG* IpeG^„ 



cqr 



(10) 



(11) 



The equilibrium GF's can be found by using the self- 



,R(A) 



hVjdTkG^ 

\R(A) 



(k,w)/(27r)™, where 



energy, 

only the imaginary part of S^q^'' should be calculated 
since the real part can be combined with the Hamilto- 
nian Hq introducing corrections to the eigenstates and 
eigenenergies of Hq that vanish as we take the strength of 
disorder to zero. Using Eq. (|7J|, we can write ImS 
^fhV'j up to the lowest order in V where 



R(A) 



7M 



N 

E 

77=1 



(12) 
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Figure 1: The SJC diagrams for atj expressed in the band 
eigenstate basis described by the indices m and n (m 7^ n). 
The bold lines correspond to G^ A ^ (that can be replaced by 
disorder-free GF's for calculation of the SJC while dashed 
lines correspond to disorder strength V. By iterating Eq. (|l(jp 
and expanding as in Eq. ©, the leading-order contributions 
in Eq. ([SJl can be expressed as two components of SJC. The 
diagrams beginning and ending with velocity involving single 
(multiple) band(s) are termed as intra(inter)-band diagrams. 



is determined solely by the electronic band structure, and 
the integral runs over the wave- vector surface correspond- 
ing to energy hx> (we will only need 7 at the Fermi level) . 

We next rewrite Eqs. © and © through GF's G^ (A) 
using Eq. ([7]). In the limit of vanishing disorder, the 
first term corresponding to vertex corrections in Eq. © 
vanishes. In the remaining terms, it is sufficient to use the 

We can identify 



diagonal GF's G^ instead of G 



eq 



the IC by combining j u with several terms from j 1 : 



2h 



d n k duo 
'2w) n 2tt 



{n F Tr[(6%v c d u G% - d u dfC e 6fj (t 



2Gfv c Gj 



Gd v c G d 



GjVcG' 



{Vc)i } 



Using integration by parts and keeping only zeroth-order 
terms in V, we arrive at Eq. ©. 

To obtain the remaining terms in Eq. (jSJ) up to the 
zeroth order in V we expand Eq. (|10p into an infinite 
series, furthermore substituting this series into Eq. (JSJ. 
In the band eigenstate representation, we can further re- 
place the GF's via diagonal ones according to Eq. ([7]). 
The resulting infinite sum has the terms of order V -1 : 



e 



N 

E 

77=1 



d n - l k Va d hkj e v - [WPjU], 

-OhkiSr, 



(27T) 



2 \dhkSr,\ [jc 



leading to the symmetric part of the conductivity, which 
describes the anisotropic magnetoresistance. The more 
interesting terms contributing to the AHE appear at ze- 
roth order in V and can be graphically represented as two 
sets of diagrams (see Fig. [T]). The inter-band diagrams, 
corresponding to calculating p = P/V + 0(V) in Eq. ([5]) 



up to the most singular (i.e., V -1 ) order, lead to Eq. ([3]). 
The intra-band diagrams, corresponding to calculating p 
in Eq. ((5]) up to the zeroth order in V, lead to Eq. Q. 
Here, P is an N X N matrix given by the solution to 
Eq. ([5]), which corresponds to the leading-order vertex 
correction to the bare velocity captured by T. 

Application to Rashba and Luttinger models. — We first 
apply Eqs. ©-(HI) to a Rashba ferromagnet with {— fl; f2} 
band gap at k = arriving at expressions (Table I in 
supplementary material [ill]) consistent with the previous 
works [8|, [9( . The vertex corrections lead to important 
contributions and should in general be considered. 

However, for inversion-symmetric systems with 
-ffo(k) = Hq(— k), the vertex corrections vanish for short- 
ranged disorder as can be seen by inspecting the P inde- 
pendent term in Eq. ([5]) . Similar vanishing of the vertex 
corrections takes place in calculations of the anisotropic 
magnetoresistance and SHE [l^ j. We apply our theory 
to 4- and 6-band Luttinger (inversion-symmetric with 
P = 0) Hamiltonians with anisotropic Luttinger param- 
eters relevant to III-V semiconductor compounds. The 
spherical model Hamiltonian in the presence of splitting 
due to interactions with polarized Mn moments can be 
written as follows within the mean-field description [l3^ : 



Ho 



h 2 

2m e 



Ti 



72 ^- 2 72 (k-j) 



-nm s , (13) 



where j is the angular momentum operator for J = 3/2, 
s is the spin operator which has to be projected onto 
the J = 3/2 total angular momentum subspace (s = j/3) 
within the 4-band model, 71 and 72 are Luttinger param- 
eters defining the light- and heavy-hole bands with the 
effective masses rriih/hh — m e/ (71 ± 272), in terms of the 
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Figure 2: SJC and IC to the AHE as a function of the 
mean-field; (a) A so — ► oo corresponds to 4-band model and 
A so = 341 meV corresponds to GaAs host, the hole density 
is 0.35 nm~ 3 , (b) the plots correspond to the In/Ga,As hosts. 
The hole densities are 0.1 nm~ 3 and 0.35 nm~ a for the for- 
mer/latter. By arrows, we mark the saturation mean-fields 
n = 25/122 meV for the In/Ga,As hosts HQ. 



free-electron mass m e , m is the magnetization polariza- 
tion direction and f2 is the mean field proportional to the 
average of local moments. For fully-polarized Mn spins, 
m is uniform and Q = NMnSJpd, where Nmu is the den- 
sity of Mn ions with spin S = 5/2, and J p d = 50 meV nm 
is the strength of the exchange coupling between the lo- 
cal moments and the valence- band electrons (lij. The 
corresponding 6-band Hamiltonian can be found in [13 1. 
As the vertex corrections vanish, all terms involving P in 
Eqs. (jHJ and Q vanish, leading, up to linear order in 0, 
to the following analytical result for Hamiltonian (fT3"|) : 



a 5p(l- Vp) + 3(1- p 5 / 2 ) 

io {i- p )(i+p-Vp) 



(14) 



where cto = (f2e 2 /37r 2 /i 2 )-\/mhh/2fiwF and p = mih/rrihh- 
SJC is in the range from 0.3ero to oq increasing as p — > 1. 

In Fig. [5^,, we present results of our calculations for the 
spherical 4- and 6-band Hamiltonians. The parameters 
are chosen to match GaAs effective masses mhh = ni e /2, 
p = 0.16 and the SO gap A so = 341 meV. The SJC does 
not change much as we vary A so . The SJC can become 
dominant for the smaller SO gaps since the IC sharply 
diminishes eventually changing sign. 

To have a more accurate description of the valence 
bands in III-V semiconductor compounds, one has to 
introduce the third phenomenological Luttinger param- 
eter 73. This leads to band warping which has strong 



effect on the IC Our calculations show that the 

SJC is accelerated by the presence of band warping. 
In Fig. [2}}, we present results of our calculations for 
(In/Ga,Mn)As for which A so = 430/341 meV, and 
(71,72,73) = (19.67, 8.37, 9.29)/(6.85, 2.1, 2.9) 0. In 
both cases, the AHE is dominated by the SJC. We 
use densities Nmu = 0.23/1.1 nm -3 for the In/Ga,As 
host leading to saturation values of the effective field 
Q = 25/122 meV 0, E^. Taking the hole density as in 
Ref. [H (0.1/0.35 nm -3 for In/Ga,As host), we arrive at 
the AHE a yx = 16/85 fT 1 cm" 1 for (In/Ga,Mn)As. Our 
results for the IC agree with the previous calculations Q 
while the total AHE overestimates the experimental val- 
ues 14, lH which is expected as the experiments are only 



on the border of the metallic regime. 

Summary. — We formulated a theory of the AHE for 
metallic noninteracting multi-band systems with the fi- 
nal result for the IC and SJC being expressed through 
the material's electronic band structure. Our derivation 
relies on the minimal coupling with the electromagnetic 
field in Hamiltonian ([T]), which is justified when a suffi- 
cient number of bands is considered. (E.g., the side-jump 
scattering in conduction bands due to spin-orbit coupling 
associated with impurities [llj can be described within 
our approach by resorting to the 8-band Kane model.) 
In contrast to the theory of the intrinsic AHE, the elec- 
tron (hole) motion in a particular band cannot be de- 
fined solely in terms of the Bloch wave functions of the 
same band in the presence of disorder-induced band mix- 
ing. The SJC does not depend on the disorder strength 
but will generally change with the type of disorder (e.g., 
short-range impurities vs. phonon scattering) . The asso- 
ciated scattering regime crossovers can be accompanied 
by a sign change of the AHE as the IC and SJC can be of 
opposite sign. The AHE sign change has been observed in 
Fe and (Ga,Mn)As [Til]. Ac measurements, furthermore, 
can quench the SJC in clean samples at intermediate fre- 
quencies t^ 1 < uj < A. We demonstrated our theory on 
electronic band structures of the 2D Rashba and 3D Lut- 
tinger Hamiltonians. Within our simple model, the AHE 
in the metallic (In/Ga,Mn)As magnetic semiconductors 
at low temperatures is dominated by the SJC. The pro- 
posed theory can be further used in ah initio calculations 
of the AHE in wide range of available metallic materials. 
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